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SUMMARY 


і method 38 presented which allows the use of nonlinear 
section lift data in the calculation of the spanwise lift distribu- 
tion of unswept wings with flaps or ailerons. This method 
is based upon lifting-line theory and is an extension to the 
method described in NACA Rep. 865. The mathematical 
treatment of the discontinuity in absolute angle of attack at the 
end of the flap or aileron involves the use of a correction factor 
which accounts for the inability of a limited trigonometric 


series to represent adequately the spanwise lift distribution. ` 


„1 treatment of the apparent discontinuity in maximum section 
lift coefficient ig also described. In order to minimize the 
computing time and to illustrate the procedures inrolred, 
simplified computing forms containing detailed examples are 
given for both symmetrical and asymmetrical lift distributions. 
A few comparisons of calculated characteristics with those 
abtained experimentally are also presented. 


INTRODUCTION 


Unswept-wing characteristics calculated by the method of 
reference 1, in whieh nonlinear section lift data are used, have 
been found to agree much closer with experimental data in 
the region of maximum lift coefficient than those calculated 
by methods in which linear section lift curves are used. It 
appears feasible that the similar use of nonlinear section data 
would yield improved results for unswept wings with flaps or 
ailerons. The deflection of a partial-span flap or aileron, 
however, causes discontinuities in the spanwise distribution 
of the absolute angle of attack. If such discontinuities exist, 
an excessively lerge number of spanwise stations must be 
eonsidered in order to obtain a solution by the method of 
reference 1 or by any other method in which the lift dis- 
tribution is approximated by a trigonometric series. A 
different method of treatment of the discontinuity is there- 
fore desirable. Developed herein is a new method of treat- 
ment involving the use of a correction factor which accounts 
for the inability of a limited trigonometric series to represent 
adequately the spanwise lift distribution of a wing with 
partial-span flaps or ailerons deflected. This correction 
factor is obtained with the aid of reference 2 and is used in 
conjunction with the system of multipliers of reference 1 to 
obtain the induced angle of attack from the spanwise lift 
distribution. The subsequent ealculation of the lift dis- 
tribution by means of successive 5 approximations is similar to 
that of reference 1. 


The mathematical treatment of the discontinuity in the 
spanwise distribution of absolute angle of attack is only part 
of the problem involved in calculations for wings with flaps 
or ailerons. The direct use of two-dimensional lift data . . 


would indicate a discontinuity in the spanwise distribution of 


maximum section lift coefficient. Obviously, the flow about 
the wing sections near the end of a flap or aileron is not two- 
dimensional. A rational method of obtaining three-dimensional 
section data from the two-dimensional data has therefore 
been devised and is presented herein. This method is... 
substantiated by experimental pressure distributions. 

In addition to the presentation of the method of making the 
calculations, simplified computing forms are given and their 
use is illustrated by detailed examples for both symmetrical 
and asymmetrical distributions. A few comparisons of cal- 
culated results with experimental data are also given. 

This report and reference 1 are intended to supplement 
each other. The reader is therefore expected to be reason- 
ably familiar with reference 1. 


SYMBOLS 


As used herein, the term section designates the аа un 
istic of a section in three-dimensional flow. 


A ‘aspect ratio | 

EP coefficients of trigonometric series for lift’ TM 
bution l 

Cp, induced drag coefficient EN 

Cr wing lift coefficient 

С, rolling-moment coefficient | _ 

Ce induced yawing-moment coefficient κ. 

Е edge-velocity factor for symmetrical part of lift — 
distribution 

E' edge-velocity factor for antisymmetrical part of. 


Iift distribution 
F factor used in altering tw чаш lift curves . 
R Reynolds number 
S pressure coefficient 
Total pressure— Local static pressure 
Dynamic pressure ) 


а, section lift-curve slope per degree 
b span of wing 

local chord of wing _ 
C, root chord 
с mean geometric chord (5/4) 
Ca, section induced-drag coefficient ке 
δι section lift coefficient 


· Ворегведсв Ела D 2283, «Method for Calculsting L Lift Distributtons for Unswept Wings with Flaps or A ferons by Use of Nonlinear Section Lift Data” by James C. Sivelisand 


Gertrude C. Westrick, 1 


628 


REPORT 1090 —NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


stall margin 

maximum section lift coefficient | 

maximum two-dimensional lift coefficient 

increment in maximum lift coefficient due to flap 
deflection (two-dimensional data) 

section lift coeficient at either end of flap or 
aileron 

section lift coefficient for part of lift distribution 
involving no discontinuity in angle of attack 

section lift coefficient for part of lift distribution 
due to discontinuity in angle of attack 

section lift coefficient for additional lift 
distribution 

section normal-force coefficient 

wing-tip helix angle generated by rolling wing 

intervals used in integration 

airfoil thickness-chord ratio . 


spanwise coordinate 

spanwise coordinate at either end of flap or 
aileron 

angle of attack, degrees | | I 

correction for induced angle of attack, degr ees 

effective angle of attack, degrees 

correction for effective angle of attack, derrew: | 

induced angle of attack, degrees 

induced angle of attack for part of lift distribution 
involving no discontmuity in angle of attack, 
degrees 

induced angle of attack containing discontinuity, 
degrees 

angle of attack for zero lift, degrees. 

angle of attack for two-dimensional lift curves, 
degrees 

angle of attack of root section, degrees 

uncorrected induced angle of attack, degrees 

multiplier for induced. angle of attack for asym- 
metrical distributions 

multiplier for induced angle of attack for antisym- 
metrical distributions 

magnitude of discontinuity, in absolute and 
induced angles of attack, degrees 

angle of twist, negative if washout, degrees 

angle of twist at wing tip, degrees 

average angle of twist 

faired angle of twist due to flap deflection 

average angle of twist due to flap deflection 

ratio of actual two-dimensional lift-curve slope to 
theoretical value of «*/90 

area multiplier for asymmetrical distributions 

area multiplier for symmetrical distributions 


„> a = s-—— e. тез — 


Tip chord 

Root chord 

multiplier for induced angle of attack for Ru 
metrical distributions 

interpolation multiplier 


taper ratio 


Um | moment multiplier for asymmetrical distributions 

T ma moment multiplier for antisymmetrical distribu- 
tions 

κ᾿ used as superscript to denote value at end of flap 


or aileron 


DEVELOPMENT OF METHOD 
= LIFT DISTRIBUTION 


The method involving the use of multipliers to obtain the 
induced angle of attack from the spanwise lift distribution 
was presented in reference 1. In this development, the lift 
distribution was approximated by a finite trigonometric 
series 


T—I 
Et! Gs Aann (1) 
л] I r 


From the 


values of ¢,c/b for each value of m, the induced angle of 
attack could then.be obtained at the points kr/r by the 


relation 
=> (F) в (2) 


where k=1, 2, 3, . r—l and Bm, denotes the multipliers 
for asymmetrical distributions. The corresponding mul- 
tipliers, X44 and ym, are used for symmetrical and antisym- 
metrical distributions, respectively. ‘These multipliers are 
tabulated in reference 1 for 7=20. The method can be used 
directly so long as there is no discontinuity in the spanwise 
distribution of absolute angle of attack caused by the deflec- 
tion of a partial-span flap or aileron. 

Lifting-line theory requires that a discontinuity in the 
distribution of absolute angle of attuck must be accompanied 
by an identical discontinuity in the distribution of induced 
angle of attack in order to avoid & discontinuity in the span- 
wise lift distribution. An analytic expression for the lift distri- 
bution associated with a discontinuity in induced angle of 
attack is presented in reference 2. The complete lift distribu- 
tion can thus be expressed as the sum of two distributions 


εις CC Cist 
To n (55) (3) 


where cı c/bô is the distribution due to a unit discontinuity 
in the induced angle of attack and c,c/b is the remainder of 
the lift distribution. These distributions are illustrated in 
figure 1. Since no discontinuity is associated with the dis- 
tribution ¢,¢/6, the corresponding induced angle of attack 
can be obtained by means of the multipliers 


au = Σ Tes ) Bmx | ч) 


where cos ΜΕ and m=1, 2, 3, . т—1. 


By definition, over the flap span 
αι = | (Ба) 
and over the unflapped span 
The total induced angle of attack is | | 
оң = оң + αι . (6) 
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If, however, the multipliers were used with the total distri- 
bution ; 


Cut 
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and Gig Were added to both sides of this equation, the result 
would be 
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Rearranging this expression gives 





r-l fec σι r—I Ci, Е 
о 2 231 вна 54 bo ) Bm | (7) 
A comparison of equation (7), for a wing with a flap or aileron, 
with equation (2), for a wing without flap or aileron, shows 
the addition of a term which is proportional to the magnitude 
of the discontinuity and acts as а correction faetor to account 
for the inability of a limited trigonometric series to represent 
adequately the spanwise lift distribution of а wing with 
partial-span flaps or ailerons deflected. For simplicity, 
equation (7) may be written as 


TE 

a, — eu, + (=) (8) 

where the uncorrected induced angle of attack is expressed as 
__ r—1 cuc 

a = 25 ( b ) в. т (9) 


and the correction factor per unit discontinuity is given by 
the equation 


miu. S fue 
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FIGURE L.— Typical lift distributions for a wing with 0.5-epan flaps. 


"The distribution ¢,¢/66 given in reference 2 may be | 
expressed in the form | 
сб 1 _ 1 — cos (6+ 8*) rent) 
= | (соз 0—сов 05) loge Tesa mt πο | (D 


α 7 
for —2=1, 0°<#<0*; for =3=0, 0+ <0<180°; and 0* and ὃ 


are in degrees. This distribution is dependent on only the 
spanwise position of the discontinuity and is independent of 
aspect ratio and taper ratio. The correction factor per unit 
discontinuity a4/5 given by equation (10) is therefore & 
function of only the spanwise position of the discontinuity. 
The distribution of =;,,/8, for which equation (11) applies, is 
illustrated in figure 2 (a). (Fig. 2 shows the method for 
obtaining various distributions of c, /6 which will be discussed 
later.) Values of σιε[ῦδ corresponding to the type of dis- 
tribution illustrated in figure 2 (а) are presented in table I 
for various values of 2y/b and 2y*/b. These values. are  . 
plotted against 2y/6 in figure 3 for even increments in 2y*/b ` 
and against 2y*/b in figure 4 for even increments in @=cos™ 2. 
In table П are given values of αε(δ corresponding to the 
values of εισ/δδ of table I. 


The distribution of αι/δ illustrated in figure 2 (a) would 
be applicable only for a wing with one outboard flap or Ἢ 


-ἰ ` O [ 
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(a) Asymmetrical distribution. 
(b) Symmetrical distribution. 
(c) Antisymmetrical distribution. 


FIGURE 2.— Method of combining distributions of a; fë. 
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aileron deflected. For a wing with symmetrical inboard 
flaps, the distribution would be obtaimed as illustrated in 
figure 2 (b) and the values of ας/δ would be obtained in like 
manner. For example, for a wing with flaps extending from 
94,* * 

—1-—0.6 to 0.06, the values οἵ æ./ô to be used would 


* 
be obtained by subtracting the values of αι/δ for Ἐν 0.6 


* 
from those for 2. ==— 0.6. Since the resulting distribution 


would be symmetrical, the values of «,/§ would also be sym- 
metrical. Similarly, an antisymmetrical distribution of 
e, fŠ and the corresponding values of a./6 would be obtained 


as illustrated in figure 2 (c). 
It should be noted that two values of а„/8 exist at the end 


* 
of the flap, one corresponding to "V о (flap side of 
* 
29*í[b) and one corresponding to “Lo (unflapped side). 


* 
Only the values corresponding to m are given in table 
II. For the unflapped side, 


| (s -(& 
б и б 


These values have no practical significance unless the dis- 
2y kx At 


continuity occurs at one of the stations "р 98 


— 1 
# 
3 0 


these stations, either of the two values may be used so long 
as the value is used with the proper section lift curve. 

In the calculation of the induced angle of &ttack by equa- 
tion (8), the values of а,,/5 must be multiplied by the magni- 
tude of the discontinuity 6. The value of 6 to be used is 
obtained from the section lift curves at the lift coefficient σι”. 
If the discontinuity occurs at one of the spanwise stations 
4) 
=-= cos e the value of (c;c/b)* is obtained as one of the 
values of c;c/b computed. If the discontinuity occurs at 
some other position, the value of (¢;c/6) * must be interpolated. 
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ΤΊ“: Ἀν 3.—Lift distribution per degree of discontinuity ір induced angle of attack for various 
spanwise locations of the discontinulty 2y*/b. - 


This interpolation may be accomplished in the following 
manner. Even though the spanwise lift distribution is con- 
tinuous through this point, the point is singular. Its singu- 
larity is due to the singularity of the corresponding point 
of the distribution of c;,c/bó which has an infinite slope but 
zero radius of curvature. At this point, from equation (11), 


Cue : x 0*sin θ᾽ | 
iJ = 8100 (12) 


Since both the distribution of e;e/b and the distribution of 
ei e/b6 contain the same type of singularity, a curve with εισ/ὂ 
as the dependent variable and c;c/bó as the independent 
variable would have no singular point and could be approxi- 
mated by a polynomial for which Lagrange’s interpolation 
formula (reference 3) is applicable. From Lagrange's formula 


(55) -z (5) >. NT 
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where 





(14) 





The number (n+1) of terms retained determines the degree - 


(n) of the polynomial used in the approximation. It has 
been found that using four terms (two on each side of the 
desired point), which define a third-degree curve, gives 
values at the end of the flap which agree within about | 
percent with the values obtained by means of the method 
of reference 4 but with about 45 terms of a trigonometric 
series. Values of the interpolation multipliers determined 
in this manner are presented in table IH for various spanwise 
positions of the end of a flap. This type of interpolation 
is suitable only if the variables are single-valued in the range 
of the interpolation. This limitation precludes the determi- 


2y* š "NEST 
nation of multipliers for =у—=0.1 but this position is out of 
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FIGORE 4.—Values of the lift distribution per degree of discontinulty μι induced angle af 
&ttack аз а function of the spanwise location of the discontinuity 10. 
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The method of determining the lift distribution by means 
of successive approximations is the same as that of reference 
| with the added step of determining the correction factor 
αι For a given geometric angle of attack, the lift distribu- 
tion i$ assumed. The uncorrected induced angle of attack 
is determined by equation (9) for asymmetrical distributions 
or by а corresponding equation using àma for symmetrical 
distributions. The multipliers уь for antisymmetrical dis- 
tributions can be used only if the lift curves are linear and 
the value of 6 is independent of lift coefficient. The value 
of (c;c[b)* is obtained by means of equation (13) and is 
divided by (c/b)* to determine с;*. At this value of c;*, ὃ is 
read from the section lift curves as the difference between 
the effective angle of attack for the section without the flap 
and that for the section with the flap, both at the station at 
the end of the flap. This value of š is multiplied by the 
appropriate values of o,/8 from table IT to obtain a. Тһе 
values of о, are added to the uncorrected values a, to obtain 
the values of induced angle of attack which are then sub- 
tracted from the geometric angle of attack to give the effee- 
tive angle of attack at each spanwise station. At each value 
of effective angle of attack, the corresponding value of c, is 
read from the appropriate section lift curve and multiplied 
by the value of οὐ at that station to obtain a check value of 


ec/b. If the check values do not agree with those originally © 


assumed, & second approximation is made and the process 
is repeated. Further &pproximations are made until one is 
found which is in agreement with the check values. 


DETERMINATION OF THREE-DIMENSIONAL SECTION DATA 


If the two-dimensional lift data were used directly, а 
discontinuity in the spanwise distribution of maximum lift 
coefficient would be indicated at the end of the deflected 
fiap or aileron. Inasmuch as this discontinuity does not 
exist in three-dimensional flow, some means of altering the 
two-dimensional data must be used to obtain what may be 
ealled three-dimensional section deta. One method of alter- 
ation is described in reference 5 for use with linear section 
lift data. In order to use this method, the lift distribution 
must be broken up into the basic Lift distribution due to 
twist and the additional lift distribution due to angle of 
attack. When nonlinear section lift data are used, the lift 
distribution cannot be broken up in this manner and the 
method of reference 5 is not applicable. For this reason, the 


TABLE HI.—INTERPOLATION MULTIPLIERS >, FOR 
с ж 
OBTAINING VALUES OF (=) AT THE END OF 
THE FLAP 


0 
0 
G 
6 
0 





method hereinefter described has been developed. Although 


— 


this development is not rigorous, the reasoning behind it із _ 


substantiated by experimental pressure distributions. 
The foundation of the method used herein is the calcula- 


tion of two curves for which linear section lift data are used. - 
One curve, the spanwise distribution of section lift coefficient 


for & wing with flaps but no other twist, is calculated by the 
method described herein but adapted for linear lift curves in 
а manner similar to that of reference 1. The other curve is 
the additional lift-coefficient distribution (constant. absolute 
angle-of-attack distribution) calculated as in reference 1. 


These curves may be calculated for any convenient angle of - 


attack and for any convenient value for the discontinuity in 
angle of attack at the end of the flap. A typical set of curves 
is shown in figure 5. It is reasoned that the root section of a 
wing with partial-span flaps would be acting most nearly like 
that of a wing with full-span flaps so that the additional lift- 
coefficient distribution is multiplied by & suitable constant 
to give the same value at the root. Likewise, it is reasoned 
that the tip section would be acting most nearly like that of & 
wing without flaps so that the additional lift-coefficient dis- 
tribution is multiplied by another constant to give the same 
value at the outermost station used in the computations (the 
0.9877-semispan station when 10 points on the semispan are 
used). The differences between the values on the curves 
(Le., ¢;(8=10°)—O.713e,,, inboard of flap end, and 


€;(6— 10^) —0.086c;,,, outboard of flap end) are then divided | 


- by the difference between the values on the additional 


lift-coefficient curves at the end of the flap (0.713¢,,— 
0.086c,,,) to obtain the factor F shown in figure 6 for several 
wings. In order to interpolate for the points on the additional 
lift-coefficient curves at the end of the flap, the multipliers 
μα, presented in table IV, have been determined. These mul- 
tipliers were obtained by using @ as the independent variable 
in equation (14) instead of c,¢/bé. It can be readily seen that 
the factor F is independent of the angle of attack and 
magnitude of the discontinuity used in the original computa- 
tions. The factor F is also relatively independent of lift- 
curve slope so that it can be used near maximum lift where 
the lift curves are nonlinear. 
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Ftacgx 5.—Typical curves used to obtain factor for altering two-dimensional data. 
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The maximum lift-coefficient values are altered by means 
of the factor F according to the relation 


δικα. m (Etnas) HP (Bei uaa) (15) 


The values of c; and o are then altered according to the 
equations νυ. | 





(Οι}οτεετεά n Cna _. -= (16) 
(Ci)unattered (61 ες) ο 
4 Eg _ Ci naz Я — (17) 
ой Qt. (οι n 


The edge-velocity factor E is used in the same manner as in 
reference 1. The value of E given in reference 6 18, however, 
probably more accurate than the ratio of semiperimeter to 
span used in reference 1. From reference 6 for unswept 


wings, u 
4 


The foregoing description of the alteration to the two- 
dimensional section data has been limited to wings with 
symmetrical partial-span flaps. For wings with deflected 
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Fig URE 6.— Factor for altering two-dimensional date for several wings. 


TABLE IV.—INTERPOLATION MULTIPLIERS jpn FOR 
OBTAINING VALUES OF 12 AT THE SPAN WISE 
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ailerons, a similar method based upon the same reasonin | 


could readily be devised. 
The reasoning behind the mothod of altering the two- 


dimensional data is substantiated by the results shown in 
figures 7 and 8 and obtained from experimental pressure 
distributions for a wing of aspect ratio of 6.0, taper ratio of 
0.5, and NACA 64-210 airfoil sections.  Chordwiso pressure 
distributions were obtained at six spanwise stations. Two 
spans of flaps, 0.496 and 0.516, were tested so that the pres- 


sure distribution at ο was obtained, in one caso, just 
outboard of the end of the flap and, in the other case, just 
inboard of the end of the flap. As shown in figure 7, the 
pressure distributions at each station are very similar for 
the two flap spans, even at the station 0.5b/2. Furthermore, 
there is а gradual spanwise change from the type of loading 
associated with an sirfoil section with flaps to the type 


‘associated with a plain airfoil. 


For comparative purposes, the section normal-force coof- 
ficients, obtained by integration of the chordwise pressure 


2y/b-09585 7 5 » J 
















v x/c*O 5.8.71 
o x/c*O S5*8.95 


v x/c«OÓ $=12.63 
ο x/c*O 5-ι2.90 
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FIGURE 7.—Experimental pressure distributions at an angle of attack of 13.1° for в wing having 
NACA 64-210 airfoil sections, an aspect ratio of 6.0, a taper ratio of 0.5, and partil-span 
split flaps. 
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distributions at each spanwise station, together with cal- 
eulated values of section lift coefficient, interpolated for the 
same spanwise station, are shown in figure 8. The calcula- 
tions for the flaps-neutral and full-span-flaps configuration 
were made according to reference 1 and those for the partial 
span-flaps configuration were made as described herein. For 
the experimental data, jet-boundary and stream-angle cor- 
rections have been applied to the angle of attack, but no 
corrections have been applied to the values of c, for the 
effects of the model supports and the boom containing the 
pressure tubes. Although some of the disagreement shown 
may be attributed to the assumptions involved in the calcu- 
Jations, most of it is believed to be due to experimental inac- 
euracies both in the two-dimensional data used in the calcu- 
lations and in the three-dimensional data shown herein, 
inasmuch as similar disagreement is evident for the full-span 
configuration as for the partial-span configuration. In 
general, however, the trends indicated in figure 8 serve to 
substantiate the method used to alter two-dimensional data. 
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WING CHARACTERISTICS 


Numerical integration.—In the numerical integration de- 


scribed in reference 1, the multipliers nw, зш, σα, ànd One 


were determined by harmonic analysis. For example, it 


was shown that 


μι ο ολα, 


be seen that equation (19) could also be obtained through 
the application of the trapezoidal rule (reference 3) where 
F(2y/6),, sin (mer) is the function of mx/r and the interval 
between points is r/r. Although equation (19) may be 
expected to give good results when f(2y/b) is approximately 


(19) .. 


-—- 


-oae а 


= ——W 


— ..— 


- oe 


т. MT ; Я M καν 
where 5, sin. —— was designated as т. By inspection, it may 
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MEL EM 


elliptic, in general the application of Simpson's (parabolic) ~~ 


rule to this problem has been found to give better results _ 


for a wider variety of curve shapes. Since г was originally ^ - 


assumed to be even, the interval m—0 to r may be divided 
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(b) 0.49-єрап split flaps. 
(d) Full-span split flaps. 


FiGC&RE ¥.—Comparizon of experimental and caleulate-t section eharateristtes tora wing having NACA 64-210 airfoil sections, an aspect ratio of 6.0, and а taper ratio of 0.5. 
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into r/2 regions, the area of each of which is 
т 241. ἐπ ὃ) (¿+ 1)r 
zb sin T Eug (ο) п T 


For the entire interval m=0 to r, 


ИЄ: (=F J S) gB- 1)”] sin 77. 


{+1 


m=1 
r—1 9 
E y). = 
where 7 is defined and used herein as 
T —f— 1\%1 ет ma a Я 
r= [8—(- D^] sin = (22) 


The relations between nm, nme, Om, 2nd c4, are the same as іп 
reference 1; namely, i 


05:772 (m #5) (232) 
Nm 1m (m=) (23b) 
гъ == cos ы (24) 
σπιαξ-2σπι 

=n, COS mn - (25) 


Values Of nms 7154, Cm, 8nd oma are given in table V for r—20. 
TABLE V.—WING-COEFFICIENT MULTIPLIERS 


9 
8 
7 
6 
δ. 
4 
8 
à. 
1 





For cases for which the end of the flap or aileron is at a 
point where m is not an even integer, а special form of 
equation (20) is required which gives additional multipliers 
for numerical integration in the vicinity of the end of the 
flap or aileron. Such a case is illustrated in figure 9 and 
the individual areas are | 


T at). 








s(2a-]-3 
Area l=" enfe Ра PEE (26) 
S T 8 (t-- 38) Gcr ((21-- 38) 
Aree 277 | παρ ge er eque Л | ση 





FiauRE 9.— Definition of symbols used [n numerical Integration for unequal increments In 


‚ Independent variable. 
5 = 
Ares 1 κ ο j, (E10? p, ο η 
(28) 
‘where s and £ are expressed as fractions of w/r. It can be 


seen that equation (28) is the same as expression (20) when 
g=t=1. Ап example of the computation of these special 
multipliers is given in table VI for 0.5-span flaps. For the 
area B in the sketch accompanying this table, equation (28) 


is used with s=1 and => and for the area O, equation (26) 


is used with js and £—1. Values of the $4, multipliers 


3 
are given in table VII for various locations of the end of a flap. 
Values of the multipliers лм, Cm, and Gne can be readily 
obtained from the values of nm, by means of equations (23) 
to (25). . : 

It should be noted that two multipliers are given at the end 
of the flap. For distributions which are continuous through 
the end of the flap, such as the distribution of c;c/b, the sum 
of the two multipliers may be used. For distributions which 
ате discontinuous at this point, such as the distribution of 
Ca C] b, each multiplier should be used ο. with its 


y* 
appropriate value. For example, for ^ 2Y 0.5 .5, the multi- 


plier —0.03023 should be used with the value at the outboard 
end of the flap, and the multiplier 0.03930 should be used 


with the value at the inboard end of the unflapped span. 


Wing coefücients.— The formulas derived in reference 1 
are repeated herein for convenience. For asymmetrical 


distributions 
| =A 5 (2 (29а) 
For symmetrical distributions - 
C A δὴ (s) m (29b) 
For asymmetrical distributions 
Co = j cue. > 2: αι). m (808) 
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For symmetrical distributions 


πα 72 c Εις ` 
For asymmetrical distributions 
(24S SEN 8 (318) 
maj 
For antisymmetrical distributions 
Στ! етс 
С.=—А >! Ἔ σπα (81b) 
т | m= 
For asymmetrical distributions 
CC 


The values of 1m, ἥπε, σα, 200 сас should be used as defined 
herein, and the summations should also include the special 
multipliers in the vicinity of the end of the flap or atleron. 


ILLUSTRATIVE EXAMPLES 

SYMMETRICAL DISTRIBUTIONS 
The method described is applied herein to a wing, the 
geometric characteristics of which are given in table VIII. 
The section lift curves are shown in figure 10 as dashed lines 
before being altered and as solid lines after being altered. 
The unaltered data are obtained by interpolation of two- 
dimensional data for the proper Reynolds number and airfoil 
thickness-chord ratio for each spanwise station їп а manner 
similar to that of reference 1. The altered curves are ob- 
tained as described previously and are used in the calculation 
of the lift distributions in tables IN and X. These tables 
were designed for use with caleulating machines capable of 


performing accumulative multiplication. The mechanics of 
computing are explained in the tables, but a few items need 


additional explanation. 

The initial approximation of ¢,c/6 is obtained in table LX. 
The computations in this table are based on the method of 
reference 6, the pertinent equations of which are modified 
to suit the present purpose. The spanwise stations, at 
which values of c;c/b are computed, are listed in column (1). 
Columns (2) to (4) are used to obtain the additional lift 


distribution e;,efe for C;,—1 according to the approximate 


relation 


(33) 


TABLE VI.—CALCULATION OF AREA MULTIPLIERS WITH 
SPECIAL CONSIDERATION AT THE END OF A.0.5-SPAN FLAP 


" Factors for areas tndleated 
I БЫ тлы u 
1 1 
0 10 s s Q. 1357080 0. 05236 
τα & 
158 9 s s . 155116 . 20686 
I | 10 18 
- 309 8 3 27 3 .119392 . 10513 
u 32 34 
454 7 = 7 | . 199959 . 16588 
| ajz H 
. 900 ὃς 9 15 15 . 136035 . 00907 
IL 1 20 
.588 6 = x 4 . 127080 ‚©0418 
4 4 178 
. 107 5 Tas £ 135 -111072 ‚1481 
1 1 2 
. 809 & 3 i $ . 092329 . 06155 
4 4 
„891 3 = = „01313 . 09508 
š 3 
1 1 2 
. 951 2 EE £ „048540 08236 
es f | ` : Ξ „024573 | _.032;6 





ο 156 309 454 588 707 8609 951 


891 
2 y/b | 


2. 
TABLE VII.—AREA MULTIPLIERS я, FOR VARIOUS SPANWISE LOCATIONS a. OF THE END OF A FLAP 
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TABLE VIII.—GEOMETRIC CHARACTERISTICS OF EXAMPLE WING 


Taper tallo, Жаан адаан ынанан Lal ieee ee ade eee” 0. 400 
Aspect ratlo, Se хасса 021 
Span (SRR πρ EGER atc a lt L... αν αν σα T. UT ҮСТ 
Root ani ph Sa a nh cm a EN E P E ане IM ao Cc da Go io ΞΗΝ arm Ma umi tim m ο 2. а iis rii 
Wing Reynolds r number, Р. „-.-..------—2-.----—--—---а-—--——----- җ----——----=-- 4.44106 


X 
a 


б. 98 
5. 43 
4. 87 
4. 35 
8. 87 
3. 44 
3.08 . 
2. 18 
2. 57 


ә P 
8 Б 


— 


For tapered wings with straight-line elements from root to construction tip: 
c ay 
a 1—(1-—А) b 


(Alter values of c/c, near tip to allow for rounding.) 


Columns (5) to (7) аге used to obtain faired values of the 
twist єз due to flap deflection according to the relations 


ay <1) (34b) 





This fairing is illustrated in figure 11. Columns (8) and (9) 
are used to obtain the average twist є; according to the 


equation 
mw 22 


The value of e; obtained by this summation is given at the 
bottom: of the table. Columns (10) to (12) are used to 
obtain an angle of attack α’ defined, for the flapped span, as 





65 


ας (at) +Hu— ita 


(36a) 
and, for the unflapped span, as 


1 


αρ τς (e;— atata 


(36b) 
The values of the geometric angle of attack а in column (11) 
include the values of continuous twist e Values of Cran in 
column (13) are obtained for the section lift curves at the 
angle a’. In the nonlinear range of the section lift curves, 
different values of οι (ш, are found corresponding to the two 
values of o/ at the end of the flaps. For the purpose of this 
table, the average value of οι τι} is used. Finally, columns 





Root [7. a» ы ыз ины ызы t a ИВОНИН ВЕРТЕР EE ORT s МАСА 01- 210 
Tipsection ана ee erimus concep une niem an μον uM nn Era Meu ώμο NAQA 61-210 
. Geometria twist, €n бөк... а-аа siluis al ως ui scm mmo m nor: ak 00 
Edge-veloci factor, Pete New eon Sea ue ο usa ou eshte ПИРЕ АТИНЕ 1. 024 
ἘΝ factor, B ————— oe Á. ланнган L 094 _ 


^ DP руа раро ро pa ро ро 
| ВЕ SZEBSSEHHS 


- c -^ ав -.- 


(14) and (15) are used to obtain the initial approximate 
values of e;c/b according to the relation 


‘ex€ _ a? 
(9) xan cy cm (37) 
which was adapted from the equation 
i ο απ _ Cial д 482 
ΖΕ τ +26 (— 9] ӨЗ) 


given in reference 6. 

The value of ë used in table IX should correspond to the 
value of οι” obtained at the bottom of the table. Some value 
of ὃ must be selected, however, for the computations needed 
to determine the value of ¢,*. It is therefore advisable to 
omit the computations in columns (10) to (15), except the 
bottom row, until а value of à has been found by trial and 
error which is consistent with the value of c;*. 

. If the lift distribution is to be obtained at a relatively high 
angle of attack, as in the example shown, the values of a 
may be in the rapidly curving part of the section lift curves 
or even greater than the angle of attack for maximum 
section lift coefficient. For the purpose of this table, the sec- 
tion lift curves may be extended in the general dircction of 
the nearly linear part of the curves in order to obtain the 
values of cnan. This procedure is justified inasmuch as 


‘the subsequent operation in column (15) will reduce the 


value of сі. 

The approximate values of cæ/b (column (15), table 1 X) 
are used m column (8). of table X and the computations 
indicated are performed. Except for the computation of ας, 
the computations in this table are the same as those in 
reference 1 for wings without flaps. The value of 6 used 
in column (16) is obtained at the section lift coefficient ¢,* 
corresponding to the value of (c;c/b) * obtained by means of 
the interpolation multipliers у». 

In this method of successive approximations, the value of 
ec/b in column (22) will usually not cheek the initial 
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FIGURE 10.—Unaltered two-dimensional and altered section Lift curves for example wing with 0.¢-span flaps. 


approximate values in column (3). The values to be used in 
subsequent approximations may be found by the equation 


A’ (+). == > ΚΑ (2) Q (39) 
Кш 

where A’(e,c/6) is the increment to be added to the approxi- 
mate values to obtain succeeding approximate values, 
A(e;c/b) is the difference between the check values and the 
approximate values (column (22) minus column (3)), and 
K and K, are constants for any partieular wing. Equation 
(39) is derived in the appendix, and values of A and K, for 
r—20 are presented in figure 12 as functions of 4Е/т. These 
values compare favorably with those empirically determined 
in reference 1 (K,=3, Ai=1, Kipi=0, апа A=8 to 10). 
Although these values were obtained for elliptie wings, they 
can be used for wings of other plen form. The number of 
terms of equation (39) needed for any particular approxima- 
tion depends upon the convergence of the approximation; 
fewer terms are needed when the differences A(e.c/b) are 
small or when positive differences nearly cancel negative 
differences. Values of K, for values of i greater than 3 are 
small enough to be considered negligible. 


Some additional explanation of the use of equation (39) 
may be desirable. For values of m from 3 to 7, equation 
(39) may be used directly if the constants Къз are neglected. 
For values of m from 8 to 10, equation (39) may be expanded 


Gp) [ee GP) e (но (8) 
KA (55). +(K,+K2A (5). LEA (22) | 


(40a) 


[к (SE) +e ($F) нда (SE) + 


A! (4) = 


(K,+K;A (z) TAA CE) | (40Ь) 
«E, iore D) екл 
OKA (55) + Kea (5) ] (40е) 


for symmetrical distributions, since a( 2) = ACE) 


ΕΠΗ. 


TABLE IX.—CALCULATION OF INITIAL APPROXIMATION OF LIFT DISTRIBUTION FOR EXAMPLE WING WITH 0.6 SPAN FLAPS 
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TABLE X,—CALCULATION OF LIFT DISTRIBUTION TOR EXAMPLE WING WITH 0,6 SPAN FLAPS 
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FIGURE 11. Fairing ofdiscontinulty due to flap deflection for Initlal approximation. ` 





AE/y 
FIGURE 12.—Coefficlents used to obtain succeeding approximations. r»20. 


For values of m of 1 and 2, equation (39) may be expanded as 


M (55) =| Ko- -KA (55) +(K.— KA (4°) + 


ка (5) ка eo] 


(41a) 


columns (16) and (17). 


κ. -ᾱ z | Ki —K3^ (ZE) +EA (9), +ка (2) + 
| . ΚΑ (99), КА (28) | 
since Б | (5) - к (ве) 


Alter convergence is obtained in table X, the values of 
с] (column (3)) and o, (column (18)) are entered in table 


(41b) 


XI and the lift and induced-drag coefficients are determined  , 


through use of the appropriate multipliers 244. 


ASYMMETRICAL DISTRIBUTIONS 


If the angle-of-attack distribution is not symmetrical, the 
asymmetrical multipliers 844 must be used in the nonlinear 
range of the section lift, curves. Ту pieal asymmetrical 
distributions are those for a rolling wing or for a wing with 
deflected ailerons. Illustrated in table XII is the case of 
a wing with flaps but without deflected ailerons, which is 
rolling at such a rate that the tip helix angle p5/2V generated 
is 0.01 radian or 0.573?. Added to the normal angle of 


attack is an increment equal to Ei (in deg), whieh is the 


equivalent Lwist of the rolling wing. In order to reduce the 
sizeof the computing form, the table of the multipliers 8,44 is 
arranged i in the form οσα. The values of 84; for a positive 
value of 2y/b are the reverse of those for the corresponding 
negative value. Instead of reversing the values of Bm, in 
table XII, the values of cfb are written in reverse order in 
column (14), Using these values with the values of В 
gives the uncorrected angle a, for the stations listed at Ше 
bottom of the table. | 
In general, the value of ë will be different for the two sides 
of the wing, so that the appropriate values must be used with 


The values of =./ë_ in column (16) 
ате those for ay" ш, 6 from table IT, whereas the values of 


αε[δ.. in column. (17) are the negative of those for у 0.6. 


TABLE XI—CALCULATION OF WING COEFFICIENTS 
FOR EXAMPLE WING 


- κ τ. TA=9021; ce 10.00] 7 





® в] ο | 
= “ T Xar 
(table X) 
(table X) @)х(09 
0.277 6.90 1.051 
1264 £A 1.251 
' H2 4.24 1.050 
‚ 213 4.40 937 
"165 6-59 1.087 
"108 1.19 — 129 
. 086 — 1 17 — 015 
„058 ә 67 . 040 
«051 9. 47 . 1268 
"029 4.92 . 148 
1158 7.30 Kir 
1150 -—4 Τὸ — TH 
Cum AX (2) X (3) 1.590 
Co e Σο 
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Another modification must be made for asymmetrical lift 
distributions since the edge-velocity factor Æ’ should be used 
for the antisymmetrical part of the distribution (see reference 
6), whereas the edge-velocity factor E has been used to alter 
the two-dimensional lift curves. From reference 6 for 


unswept wings, 
J 16 
E” = iT 


This value may be taken into account in the following 


(42) 


manner: The symmetrical part of the effective angle distri- 


bution is 


(α-- а) (a— ax) r—k 
о 


which is used directly with the altered lift curves. The 
antisymmetrical increment in the angle distribution 18 
(а оц) а (a— ai):—x 
| 2 
which must be multiplied by the ratio H/E’ in order for it 
to be used with the same lift curve. The effective angle is 
therefore 
 (a— est (a—a);-& Е ([α--αὐε--(α--αι)... 
ksi 2 ΓΕ; 2 








--(α--αι)ε-- = = [(0— αι)ε--(α-- он), а] 
== 20—04 Ао, (48) 
where | 
Ao =Ë RE (а—афк—(«— a)i (44) 


Equation (44) is computed in column (21) of table ХН. 
Other than these modifications, the computing required 
for table XII is similar to that for table X. 


DISCUSSION 


The lift characteristics of two wings without flaps and 
with 60-percent and full-span flaps have been calculated by 
the method described herein and are presented in figure 13 
together with experimental results from reference 7. One 
wing had NACA 64-210 sections and was equipped with 
split flaps. The other wing had NACA 65-210 sections and 
was equipped with split, single slotted, and double slotted 
flaps. For the split-flap conditions, the agreement between 
calculated and experimental results is quite satisfactory, 
whereas the agreement is less satisfactory for the single- and 
double-slotted-flap conditions. Since the discrepancies occur 
for both the 60-percent-span and full-span conditions, they 
are probably due to differences between the two-dimensional 
and three-dimensional section characteristics rather than due 
to the method of calculating. Some of the discrepancies in 
maximum lift coefficient may be due to the fact that the 
characteristics of these wings were extremely sensitive to 
small surface irregularities. 

The stalling characteristics of these same wing-flap com- 
binations are presented in figure 14 together with the cal- 
culated stall-margin distributions. The stall margin Ae; 1s 
the difference between the maximum section lift coefficient 
altered as described herein and the section lift coefficient at 
the maximum wing lift coefficient. The spanwise location 
of zero margin should correspond to the location of initial 
stall, and the margin at other spanwise locations is an 
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(a) NACA 64-210 sections; split flaps. 
(c) NACA 65-210 sections; single slotted 
flaps. 


(b) NACA 65-210 sections; split flaps. 
(d) МАСА. 65-210 sections; double slotted 
fisps. 


FIGURE 13.— rimental and calculated lift carves for wings with and without flaps; aspect 
ratio, 9.021; taper ratio, 0.4; washout, 2°. Experimental points designated by symbols. 
indication of the manner in which the stall spreads. In 
generel, the agreement between the experimental and 

calculated stalling characteristics is very good. 

The foregoing comparisons between calculated and experi- 
mental results were made for the вате Reynolds number. 
If possible, such comparisons should be made at the same 
Mach number also, unless the Mach number is low enough 
to have & negligible effect. Even at relatively low values of 
free-stream Mach number, adverse compressibility effects on 
maximum lift coefficient have been noted (reference 8) when 
sonie velocity is reached locally on & wing. Similar effects 
in two-dimensional flow have not as yet been thoroughly 
investigated so that calculations based on available two- 
dimensional data must be limited to subcritical Mach numbers. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 
LANGLEY FigLD, Va., November 13, 1950. 
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TABLE XII.—CALCULATION OF ASYMMETRICAL LIFT DISTRIBUTION FOR EXAMPLE 
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о. = 14.0? G = 1.25. 


e = 15.0? G, = 1.32 


FiavRE 14.—Experimenta! and calculated stalling 
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(а) NACA 64-210 sections; split flaps. 
characteristics for wings with and without flaps; aspect ratio, 9.021; taper ratio, 0.4; washout, 2°. 
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APPENDIX 
DETERMINATION OF COEFFICIENTS USED TO OBTAIN SUCCEEDING APPROXIMATIONS | 


In the method of successive approximations to determine 
the lift distribution, it is desirable to reach convergence with 
a minimum number of approximations. This desideration 
necessitates that each successive approximation be obtained 
from preceding computations in some manner. The manner 
in which these operations were performed in reference 1 and 
the coefficients used therein were determined empirically. 


It is shown hereinafter, however, that the procedure and 


similar coefficients can be determined theoretically. 

In the following derivation, A(c;ic/b) is used to designate 
the difference between the.check values and the approximate 
values, and A'(e;c/b) is used to designate the increment to be 
added to the approximate values to obtain the succeeding 
approximate values. The section lift curves are assumed to 
be linear and $, to be zero.. The operations performed 
during the first approximation (table X) may then be 
represented by the equation 


ο πο... 


If suitable increments A'(cic/b) are. chosen so that the check 
values become equal to the approximate values for the 
gecond approximation, 


сдрн реб Gon), as 


The difference between these equations is 


(2) S AA (22), шы 9) NU (2) (АЗ) 
From reference 1 " 
n т 
8x sin — 
P 
The lift-curve slope may be expressed as . 
η т? A 
a= ep сс =e oo (AB) 
where 7 is the ratio of the actual two-dimensional lift-curve 
slope to the theoretical thin-airfoil value. For elliptic wings, 
(š jh ЕД sin Ë = ` (A6) 
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"Therefore 
(A7) 


which, for constant η, is constant for all spanwise stations. 
Substituting this value into equation (АЗ) yields 


(144E Al (ЛЕ Мак у (5) =2 aE ес 
nr m= Ё № b k 


Equation (А8) represents r/2 simultaneous equations which 
may be represented in matrix form as 


(A8) 


εις _АЕГ,® σις (Α9) 


[М] [^ 


where [M] is x P with all the principal diagonal elements 
equal to 1452 and the other elements are ER This 


matrix can readily be put into a symmetrical form and its 
reciprocal obtained by one of the standard methods pre- 
sented in reference 9. Then 


" ` É ae) АЕ AE MJ- [А че 


Кога given value of AE/y, equation (А10) may be expressed 
in the form 
(5E | 
b mai 


beeause of the particular properties of this reciprocal matrix. 
For conyenience K, is made equal to unity and the values of 
К and K, are adjusted accordingly. The values given in 
table XIII and figure 12 were obtained for various values of 
AE/n and r=20. Ἢ 


(A10) 


(2) = Le КАА (А11) 


TABLE XIII -—COEFFICIENTS USED TO OBTAIN 
SUCCEEDING. APPROXIMATIONS 
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